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Scaling dimension of fidelity susceptibility in quantum phase tran- 
sitions 
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Abstract. - We analyze ground-state behaviors of fidelity susceptibility (FS) and show that the FS has its own 
distinct dimension instead of real system's dimension in general quantum phases. The scaling relation of the 
FS in quantum phase transitions (QPTs) is then established on more general grounds. Depending on whether 
the FS's dimensions of two neighboring quantum phases are the same or not, we are able to classify QPTs 
into two distinct types. For the latter type, the change in the FS's dimension is a characteristic that separates 
two phases. As a non-trivial application to the Kitaev honeycomb model, we find that the FS is proportional 
to L 2 InL in the gapless phase, while L 2 in the gapped phase. Therefore, the extra dimension of InL can be 
used as a characteristic of the gapless phase. 
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Introduction — Let us consider a general QPT [1] occurring 
in the ground state |*Po(/l)) of a quantum many-body Hamilto- 
nian 

H = H + AH It (1) 

where Hj is the driving Hamiltonian and A denotes its strength. 
The motivation of the fidelity approach [2] to QPTs is that 
the overlap between two ground states separated by a small 
amount 5A, i.e. F(A,A + 6A) = \^¥ (X)\V (A + 6A)}\, is ex- 
pected to show a minimum around the critical point A c due to 
the dramatic change in structures of the ground-state wavefunc- 
tion [3,4]. This interesting insight to QPTs from quantum infor- 
mation theory [2] had then been demonstrated in a few strongly 
correlated systems [5-8]. It was realized consequently that the 
leading term of the fidelity, called the FS [9] or the Riemannian 
metric tensor [10], should play a key role in such a new ap- 
proach to QPTs. After that, various issues based on the fidelity 
or its leading term [11-27] , including scaling and universal- 
ity class [13, 14], and its role in topological QPTs [23-26] etc, 
were raised and addressed. 

However, it seems to us that all relevant studies took it for 
granted that the FS, in general quantum phases, has dimension 
d, \.&.,Xf k L d where d(L) is the dimension (length) of the 
system. In this work, we will show instead that the FS has its 
own dimension depending on both the scaling dimension of the 
driving Hamiltonian and long-range behaviors of their correla- 
tions. The distinct dimension of the FS means that, in a class of 



quantum phases, the adiabatic response of the ground state to 
driving parameter is no longer proportional to the system size. 
This property is different from our previous understanding of 
critical phenomena from statistical quantities, such as energy 
that is extensive in a thermodynamic system. Therefore, the 
observation not only provides a new angle to understand the 
role of FS in QPTs, but also is of fundamental importance to 
apply quantum adiabatic theorem to scalable systems. 

The critical exponents of the FS and their scaling relation 
then will be proposed in a more general way. Clearly, the FS's 
dimension can be changed, or not changed in QPTs, this prop- 
erty classifies all QPTs into two different types. For a class 
of QPTs, the change in the FS's dimension is a characteristic 
that separates two phases. As a non-trivial application, we will 
show that, in the gapless phase of the Kitaev honeycomb model, 
the dimension of the FS becomes L 2 In L rather than L 2 in the 
gapped phase. Therefore, the extra dimension of In L becomes 
a characteristic of the gapless phase. 

Scaling analysis revisited — The FS is defined as the leading 
term of the fidelity^ = - lim^^o 2 \nF I (6 A) 2 . For the Hamil- 
tonian system [Eq. (1)], the ground-state FS can be evaluated 
from [9, 10] 



Xf{X) = J] 



n*0 



\(VM\H,\W))\ 2 

[E„(A) - E (A)] 2 
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where ff(A)PF B (A)) = E„|T„(/l)>. Here we would like to point 
out that, in the low-energy spectra of the Hamiltonian, only 
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those excitations with nonzero ( x ¥ n (A)\H I \ l ¥o(A)) have contri- 
bution to the FS. So if we define A as the energy gap be- 
tween the ground state and the lowest excitation with a nonzero 
(f¥ n (X)\Hi^¥ (X)), the FS satisfies the following inequalities 



XF 



«#0 



[E n (A) - E (A)] 



1 d 2 E(A) 
"2A dA 2 



A- 2 [<*Fo(^)l^l v i , o(^)> - <W)|///| W)> 2 ]. 



Let us suppose Hj = £ r V(r) and N = L d . 1) If the system 
is gapped, the FS has the same (or smaller) dependence on the 
system size as the second derivative of the ground-state energy. 
2) If Xf/N still increases with the system size, the ground state 
must be gapless. 

The above inequalities can only tell us some qualitative in- 
formation. To find the scaling relation, we resort to the defini- 
tion of the FS in terms of correlation functions [9], 



Xf_ 
L d 



rG(r, T)dT, 



where G(r, t) = (V(r, t) V(0, 0)> - (V(r, 0)>< V(0, 0)>. Under the 
scaling transformation [28] r' = sr, -f = s^t, V(r') = s~ Av V(r), 
it has been found that the FS scales like Xf/L 1 * oc i d+2 (- 2A v 
[13] where £ is the dynamic exponent and Ay is the scaling 
dimension of V(r) at the critical point, and xf is believed to 
be extensive away from the critical point. The latter "common 
sense" is also the reason that, in the most previous studied, the 
FS is believed to have the same critical exponent a on both 
sides of the critical point, that isxf /L d oc \A - A c \~ a . However, 
it is not universally true. 

For a general <f -dimensional system, the correlation function 
G(r, t) under the scaling transformation, becomes 



G(r',r') 



= G(sr,/r) 



(4) 
(5) 



where A A comes from the transformation A' = s Aa A. Let s 
r~ 1/f , we then have 



G(r,r)^T 2 ^- d ^G(rT-^,l). 



(6) 



If we rearrange the expression, 

G( rr- 1 ^, 1) = [n --i/f]2(A,-<0 /(rr -i/f). (7 ) 

which defines the scaling function f(rr~ 1 ^). The correlation 
function becomes 



1 



(8) 



where A v = d - A A is just the scaling dimension of V(r). If the 
correlation length is divergent, the correlation function decays 
algebraically, and 



Xf_ { L d+2 ^ 2Av , 2A v -2£±d 
L d ~ \ InL, 2Av -2£ = d ' 



(9) 



For some collective systems, such as the Lipkin-Meshkov- 
Glick (LMG) model [29, 30], we will show below that the scal- 
ing dimension of cr in the polarized phase is A^ 1 , hence Xf 
is intensive. Therefore, the FS has its own distinct dimension 
instead of real system's dimension. The dimension is related 
to the size-dependence of the fidelity of a given driving Hamil- 
tonian, and as we show in our another paper that L d - actually 
determines the duration time scale in the quantum adiabatic the- 
orem [31], we will call it quantum adiabatic dimension(QAD), 
and use d a to denote it hereafter. 

Therefore, to judge a quantum criticality from the FS, the 
scaling relation and conditions should be revised. If the di- 
mension otxa is cC <-> above (below) the critical point, then the 
rescaled FS, i.e. xf IL d ? , as an intensive quantity, scales like 



1 



Xf 



L d i \A-A c \ a± ' 



(10) 



in the thermodynamic limit. Since the FS usually shows a max- 
(3) imum around the critical point and scales like Xf\a=a c x L m 



(p. = 2d + 2£- 2Av) [13], the general scaling relation should be 



a = 



(ID 



Here v is the critical exponent of the correlation length, which 
defines the length scale of the system around the critical point. 
In case the exponent v is also different (though we did not find 
such a case yet), a more general relation is a* = (ji - d^)/v ± . 
Hence the condition of singularity in the FS around the critical 
point is 

H > dt (12) 

instead of fi > d [13]. A further remark is that, even if the 
equality condition of Eq. (12) is satisfied, it is still possible for 
the system to undergo a logarithmic divergence at one side of 
the critical point, i.e. 



(13) 



if XF/L da \A=A c diverges as InL. Therefore, if and only if the 
QAD (including logarithmic dependence on the system size) 
of the FS are the same above, at, and below the critical point 
A c , it is firm to say that the FS does not have singular behavior 
around the critical point. 

Eqs (10-13) represent our first main result of this work. 
These key relations define also a criteria to judge quantum 
phase transitions in perspective of the FS. They are a gener- 
alization of the results by Venuti and Zanardi [13]. In short, in 
order to study the scaling behavior of the FS, one should con- 
sider the QAD of the driving Hamiltonian in the corresponding 
quantum phases. The QAD can be generalized to a Hamilto- 
nian defined in a high-dimensional parameter space. In this 
case, different driving Hamiltonian might have different QAD. 
In the following, we will check the validity of about analysis in 
a well studied QPT of Landau's type and a topological QPT 

Examples: the LMG model and the Kitaev honeycomb 
model — We first check the above scaling relations in an exactly 
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Fig. 1 : (Color online) LEFT: The FS Xf/L 2 as a function of In L for the 
Kitaev model along the evolution line J x = J y defined on J x +J y +J z = 1 
plane. In cases of J- < 1/2, the explicit fluctuation of the FS is due 
to appearances of infinite critical lines as the system size increases 
[32, 34]. Nevertheless, the fluctuations do not affect the logarithmic 
divergence of Xf/L 2 . RIGHT: The rescaled \J- - 1/2\ 1I2 xf/L 2 as a 
function of In \J- - 1/2| for a system of L = 20001. 



state of the Kitaev honeycomb model consists of two different 
phases [32], i.e, a gapped phase with Abelian anyonic excita- 
tions and a gapless phase with non- Abelian anyonic excitations. 
The critical behavior of the fidelity in the model was previously 
addressed by two groups [24, 25]. None of them addressed the 
QAD of the FS, but took it for granted that d a = d. We define 
the phase diagram on the plane J x + J y + J, = 1, and consider a 
certain line J x = J y along which the ground state of the system 
evolves at zero temperature. In this case, J z is the only driv- 
ing parameter. The Hamiltonian can be explicitly diagonalized 
in the flux-free subspace [32, 34] and the ground-state FS for a 
system of N - 2L 2 (odd L) sites is [25] 



Xf 



-y 

16 ^ 



sin q x + sin q x 



(16) 



where q^ - 2nn/L, n = -(L - l)/2, • • ■ , (L - l)/2, and 

e q = J x cos q x + J r cos q y + J z , 
A q = J x sin q x + J y sin q y . 



(17) 



solvable model [29,30], i.e. LMG model. Its Hamiltonian reads 

Himc =-^Z(^ +rCT '>^)- h (14) 

where y denotes the anisotropicity. The prefactor 1 /N is to en- 
sure a finite energy per spin. In the thermodynamic limit, the 
ground state of the system for y + \ undergoes a second-order 
QPT at h c — 1 . If h > h c , the system is fully polarized, while 
if < h < h c it is a symmetry-broken state. The ground-state 
fidelity of the LMG model has been addressed previously by 
Kwok et al [18]. They found that^/r is proportional to the sys- 
tem size /V in the symmetry-broken phase, while is intensive in 
the fully polarized phase. This difference makes that critical ex- 
ponents to be different at both sides of the critical point. Kwok 
et al [18] obtained that d* = for h > 1, d~ = 1 for h < 1, 
fi = 4/3ath c = 1, and the correlation length exponent v = 2/3, 
then the critical exponent above h c is a + — (ft - d^)/v - 2, 
while below h c is a~ = (ji- d~)/v = 1 /2. 

The explicit differences of d* and a* in two quantum phases 
of the LMG model has already been a straightforward demon- 
stration that the FS has its own dimension. According to Eq. 
(3), the deep reason of such a difference is that en? ~ A^ 1 in 
region h > 1, while crl ~ N~ 1 ^ 2 if < h < 1, which leads to 
that the FS is intensive in the polarized phase, while extensive 
in the symmetry -broken phase. 

As a non-trivial application, we study the FS in a topological 
QPT [33] occurring in the ground state of the Kitaev honey- 
comb model [32]. The model is associated with a system of 1/2 
spins which are located at the vertices of a honeycomb lattice. 
The Hamiltonian reads: 



In the thermodynamic limit, the FS becomes 



H = -J x J] o*jvi-J y 2 



a 3 ^ - J z 



.vbonds 



E 

cbonds 



ayr k , 



(15) 



where j, k denote two ends of the corresponding bond linked to 
a vertex, and J K (k = x, y, z) are coupling constants. The ground 



%=-^L dqx L dqy 



sin q x + sin q x 



e 2 + A 2 



(18) 



In the gapped phase J z > l/2,e 2 +A 2 > 0. So there is no pole 
in the integrand of Eq. (18). Xf/L 2 is intensive and the QAD of 
the FS is d + a = 2. However, in the gapless phase, £q+A 2 has zero 
points in k space. Expanding the integration around the pole of 
the integrand, we find^ f /L 2 oc f jdk ~ lnL,where A is a 

cut-off. Therefore,;^//, 2 manifests distinct size dependence in 
the gapped and gapless phases. As a numerical demonstration, 
we show the FS Xf/L 2 as a function of InL in Fig. l(left) for 
various J z . Obviously, Xf/L 2 does not change as the system 
size increases in the gapped phase, but is proportional to In L 
in the gapless phase. At the critical point J- = 1/2, it has been 
shown that^f oc L 5 ^ 2 [25]. Therefore, according to the scaling 
analysis, the critical exponents of the FS should be different at 
both sides of J z - 1/2. At J z - (l/2) + it has been already 
obtained a + 0.5, here we find that 



Xf\J z 



JzA m 



L 2 InL 



ln|/-7 z , 



(19) 



at J z = (l/2)~ [simplified as 1/2-ln for QAD and shown in 
Fig. 1 (right)]. Moreover, differ from both the quantum Ising 
model and the LMG model, the FS has higher dimension 2+ln 
(for^-f oc L 2 InL) than the system's dimension 2 in the gapless 
phase. The extra dimension of In L has a special meaning be- 
cause it only exists in the gapless phase with non- Abelian any- 
onic excitations. Therefore, it can be used as an characteristic 
of the gapless phase in the Kitaev honeycomb model. 

From the point view of Eq. (3), the deep reason behind the 
extra dimension In L is that the bond-bond correlation function 
of the z-bond in Eq. (15) decays exponentially in the gapped 



p-3 



S.J. Guetal. 



Table 1: Critical exponents \i,v, dp, a*, for the Ising model [4], the Kitaev toric model(KTM) [26], the LMG model [18], and the Kitaev 
honeycomb model (KHM). 



Model 


11 V 


d + a a + 


d a a 


ID Ising model(/z e = 1) 


2 1 


1 1 


1 1 


KTM(i e = iln(V2 + l)) 


In 1 


1 In 


1 In 


LMG model (h c = 1) 


4/3 2/3 


2 


1 1/2 


KHM(7 C = 1/2) 


5/2 1 


2 1/2 


2+ln 1/2-ln 



phase, while algebraically, i.e., G(r, 0) ~ 1/r 4 , in the gapless 
phase [25]. Therefore, the FS in Eq. (3) becomes 

r 

The InL divergence of the FS Xf/L 1 * means that the dynamic 
exponent should be f = 1 . 

Finally, we summarize the critical exponents of the FS in the 
above two typical models, the quantum Ising model, and the 
Kitaev toric model(KTM) in Table. I. The critical exponents 
of the FS in the latter two models are referred from Ref [4] 
and [26], respectively. The table can be prolonged as one con- 
sider more and more QPTs. No matter what transitions are in- 
cluded, the QPTs can be classified into two distinct classes. For 
the first class, the FS has the same QAD in the both quantum 
phases, such as the quantum Ising model, while for the another 
class, the FS has different QAD, such the LMG model and the 
Kitaev honeycomb model. Moreover, for the Kitaev honey- 
comb model, we can see that the QAD actually plays a role of 
characteristic that separates two quantum phases. These con- 
clusions constitute our second main result. 

A brief summary and a challenge — In summary, we have 
analyzed the dimension of the FS in quantum phase transitions. 
We have shown that the QAD of the driving Hamiltonian is not 
always the same as the system's dimension. The scaling rela- 
tion of the FS in various QPTs is established on more general 
grounds. The FS might have distinct critical exponent at both 
sides of the critical point. So the QPTs can be divided into two 
classes based on the criteria if the QAD is changed or not dur- 
ing the phase transition. Our results also show that the QAD 
can be used as a characteristic that separates two topological 
phases in the Kitaev honeycomb model. Therefore, the QAD 
provides a quite distinct tool instead of the traditional order pa- 
rameter to study QPTs. 

Clearly, the QAD provides a unique classification of quan- 
tum phases. From this point of view, a challenging problem 
might be that, for a <f -dimensional system, does there exist such 
a quantum phase with a QAD of d + 1 ? 
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